Reconstruction of Gravitational Lensing Using WMAP 7- Year Data 
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Gravitational lensing by large scale structure introduces non-Gaussianity into the Cosmic Mi- 
crowave Background and imprints a new observable, which can be used as a cosmological probe. 
We apply a four-point estimator to the Wilkinson Microwave Anisotropy Probe (WMAP) 7-year 
coadded temperature maps alone to reconstruct the gravitational lensing signal. The Gaussian bias 
is simulated and subtracted, and the higher order bias is investigated. We measure a gravitational 
lensing signal with a statistical amplitude of C = 1.27 ± 0.98 using all the correlations of the W- 
and V-band Differencing Assemblies (DAs). We therefore conclude that WMAP 7-year data alone, 
can not detect lensing. 



I. INTRODUCTION 

Gravitational lensing of the Cosmic Microwave Back- 
ground (CMB) provides information on the mass distri- 
bution between the surface of last scattering and the ob- 
server, thus potentially providing information, for exam- 
ple, on dark energy and neutrino masses. In addition, 
gravitational lensing causes E'-modes to be converted into 
large angular scale i?-modes, thereby potentially contam- 
inating B-mode signature of inflationary gravitational 
waves Because lensing deflects CMB photons by ap- 
proximately 3', a pcrturbative treatment to first order is 
generally valid. An estimator for the deflection angle has 
been devised by Hu @, Q . 

The first attempt to detect lensing by Hirata et al. 
[|| used the cross-correlation between the WMAP 1-year 
data and selected luminous red galaxies (LRGs) from the 
Sloan Digital Sky Survey (SDSS). No statistically signif- 
icant signal was found. The first detection of lensing 
was performed by Smith et al. [j| who used the cross- 
correlation between the NRAO VLA Sky Survey (NVSS) 
of radio galaxies with a higher mean redshift than the 
Sloan LRGs and a fully-optimal lensing estimator on the 
statistically more powerful WMAP 3-year data. Evi- 
dence for lensing was found at the 3Aa level. Using a 
similar estimator as in [H, Hirata et al. @ obtained re- 
sults consistent with, though at slightly lower significance 
than Q, using WMAP 3-year data, LRGs and quasars 
from the SDSS data, as well as data from the NVSS. Rc- 
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cently, Smidt et al. used an estimator based upon 
the kurtosis of the CMB temperature four-point corre- 
lation function to estimate lensing from WMAP 7-year 
data only and claimed evidence for lensing at the 2a level. 
Recently, the Atacama Cosmology Telescope (ACT) col- 
laboration successfully detected gravitational lensing Q 
at the 4cr level. The South Pole Telescope (SPT) de- 
tected the effects of gravitational lensing on the angular 
power spectrum]^]. 

In this paper we present a search for gravitational lens- 
ing using the WMAP 7-year data alone and the standard 
optimal quadratic estimator 0, which differs from the 
kurtosis estimator of ■ We apply the quadratic estima- 
tor to WMAP-7 temperature maps alone for the first time 
in the hopes that our analysis might serve as a touchstone 
allowing for consistent comparison between different lens- 
ing extraction techniques. We review the notation for 
full-sky reconstruction of gravitational lensing in Section 

II. We discuss the sky-cut used in our analysis in Section 

III. Then we introduce our modified estimator in Sec- 
tion IV making use of the optimal quadratic estimator 
of 0. We introduce the WMAP 7 -year data in Section 
V, and describe the details of the calculations, including 
the noise model, and analysis in Section VI. Results of 
a null test are shown in Section VII, and we discuss the 
conclusions of our work in Section VIII. 



II. GRAVITATIONAL LENSING 

The effect of lensing on the CMB's primordial temper- 
ature T in direction n can be represented by 



T(n) = T(n + d(n)), 



(1) 
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where T is the lensed temperature and d(n) = V</>, with 
4> being the lensing potential. The two-point correlation 
function of the temperature field following is: 
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where the second term encodes the effects of lensing with 
the weighting factor f^J, given by 

fiui' = Cf T oFl'U + Cf T a Fi L i>- (3) 

Here Cj T are the unlensed temperature power spectra, 
and 



qFiw = 
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The lensing estimator is constructed from an average 
over a pair of two-point correlations 0, Q and has the 
form 
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The requirement that the estimator in Eq. 
biased and has minimal variance results in 



AY 



L(L + 1)(2L + 1) 



and 



cTT 
JILV 



2C\ ot C\? 



(6) 



(7) 



with C\ ot = Cf T +Nf T , where Cf T are the lensed power 
spectra and Nf T is the instrumental noise. In the fol- 
lowing, the summations are from / and V = to 750 
and |m| < i, \m'\ < I'. The WMAP 7-year data do not 
contain additional information at higher multipoles. 

To reduce computation time we follow (Tol | and define 
three maps for the TT estimator: 



A T (n)=Y,-^tT lrn0 Y lm (n), 
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and take the inverse Spherical Harmonic Transform 
(SHT) of A T X and A T Y to get 



Tlm = Plo I dn +1 Y£ M A T X 
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Using Eqs. (jHJl, © and (fl0|) the expression for d^Jt in 
Eq. ([5]) becomes 
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A similar procedure is followed for the efficient calcula- 
tion of Ajr in Eq. (j6|) . The resulting expression is given 
in [ll| (originally proposed in fl2l|): 
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d(cos6>) (16) 



47 r cr+ivr 



(18) 



ff ±1 (g) = g + 1 ) C i?+ ^TT ^lW. (19) 



here d' ss / (0) are Wigner d-functions. 



III. SKY CUT 



C TT 
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In order to eliminate contaminated data, regions such 
(9) as the galactic plane and bright point sources in the full- 
sky map must be removed using a mask, thereby intro- 
ducing a sky-cut. For example, in 0, the Kp2 mask 
was used to make 84.7% of the sky uncontaminated. In 
(10) 0; the more conservative KQ75 mask was used to clean 
artifacts around the galactic plane and point sources. 
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The sky-cut can be removed as a separate component 
to get a full-sky map before we process the data. One 
such technique is the "inpainting" method in which the 
estimated values of pixels in the map are substituted for 
those removed by the mask. Perotto et al. have simu- 
lated the full sky reconstruction for PLANCK Pj. The 
full-sky map recovered in this way will bias the lensing 
reconstruction slightly. 

Another method proposed by A. Benoit-Levy [l4[ 
apodizes the masked regions of the map and inpaints the 
masked regions of the map by constrained Gaussian ran- 
dom values of the unlensed temperature. In this way, 
the sky-cut-induced coupling approximately reduces to a 
unit matrix. However, for WMAP, we have to remove a 
big portion of the sky, reducing / s k y dramatically to 0.3. 
The unbiased estimator could be scaled up by a factor 
of l// s ky, but the signal-to-noise ratio would be reduced 
significantly. This means the uncertainty of the recon- 
structed signal would be larger. 

As opposed to a separate-component solution, we ob- 
tain an all-inclusive lensing reconstruction pipeline, using 
the built-in filter of the estimator to treat the data with- 
out pre-conditioning it. The optimal estimator for the 
potential based on the maximum likelihood is derived by 
Hirata fl5j . The full inverse variance (C + N)- 1 , instead 
of (Cj T + Nj / T )~ 1 , was used by @ because it is an op- 
timal filter when there are sky-cuts and inhomogeneous 
noise. The sky-cut generates artifacts in harmonic space, 
as does lensing. (C + N) _1 can be used to filter those 
modes affected by the sky-cut. However, we do not use 
this filter because the inversion of (C + N) is computa- 
tionally challenging [f| ; instead we use the estimator Eq. 
([5]) which is identical to the one of [f| , and it is an excel- 
lent approximation to the maximum likelihood estimator. 
We note that, while (Cj T + A r J T ) _1 will be suboptimal 
to a full (C + N) -1 filter, it preserves the simplicity and 
efficiency of the lensing reconstruction procedure. 
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FIG. 1: The higher order bias calculated from {CT - iv£ 0) ) - 
Ct d for all correlations of the WMAP's W- and V-band DAs. 
The simulated higher order bias from averaging 700 (to be 
discussed in Figure I10p realizations is shown in orange. For 
comparison, the simulated lensing signal is shown in green. 
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IV. THE LENSING ESTIMATOR 



following Eq. , and a harmonic mode of the Gaussian 
bias is estimated as 



For WMAP, we modify the estimator slightly to deal 
with the instruments' anisotropic temperature noise. 
The observed lensed temperature map T is given by 

T(n) = M(n) [ J dn'T(n')B(n, n') + 7V(n)] (20) 

and likewise the "observed" unlensed temperature map 
f is 

T(n) = M(n) [ J drif(n')B(n, n') + N(n)] (21) 

Here M(n) represents the mask, i?(n,n') the beam, and 
N(n) the noise. 

For a pair of maps a and f3, "TT(a x /?)" denotes the 
cross-correlation between these two temperature maps. 
A harmonic mode of the reconstruction including noise 
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Here C are the power spectra of the observed lensed 
temperature, determined from (T; TO T;' m '). As was done 
in [8( and [16j we use the same power spectra in Eq. 
([22]) and Eq. ([23"]). In order to deal with the non- 
uniform noise distribution in the WMAP data, we sym- 
metrize d^^} QX ' 9 '' as in 0, denoting the symmetrized 
cross-correlation u TT(a» (3)" between these two temper- 
ature maps, 
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and including noise, and = {N^ m Nlm) as the Gaussian 

, „. rprp.H n bias, with the superscript "TT(a»/3)" omitted. Thus we 

lm +jn lm ^5) obtain 

We refer to C£ s * = (d* LM d,LM) a s the reconstruction 
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FIG. 2: The normalized likehood of the amplitude of the 
higher order bias limited to the region 20 < L < 170, to 
the simulated lensing signal. This confirms that the higher 
order bias is consistent with zero and negligible. 



following a reasoning similar to the one near the end of 
Section II. 

The two-point correlation of the Gaussian bias esti- 
mator is essentially a four-point correlation function of 
the primordial temperature modes. It should be care- 
fully subtracted since, for a noise-dominated experiment 
such as WMAP, the Gaussian four-point bias is sev- 



eral orders of magnitude larger than the lensing power 
spectra. In [§[ phase-randomized data maps are used 
to simulate this Gaussian bias. However, this approach 
does not work for the present lensing reconstruction since 
WMAP's noise is not isotropic. Evidence for this can be 
seen from the normalization factor A^ T< ' ax ^ which is 
not equal to j\f^ TT ^ a '^ whereas they should be equal 
for isotropic noise (Io| . The normalization factor Eq. 
(|27p only contains the partial contribution coming from 
the non-isotropic noise while the Gaussian bias squared 
from Eq . (|2~5j) consists of all the correlations generated by 
the non-isotropic noise, see [T3| and [l8[ . If the phases of 
the WMAP temperature maps are randomized in order 
to remove the lensing-induccd coupling between modes, 
it will also remove the strong correlation of the noise. 
The Gaussian bias calculated in this way will be signifi- 
cantly lower than that from the standard approach [l9j |. 
So we have to perform simulations which use the simu- 
lated WMAP noise and temperature maps, rather than 
the randomized WMAP data to get the Gaussian bias 
term. 

The deflection power spectrum is 

ct = < [ d £<~»] * - [<T" ] ] * Kit^y 

n (28) 

This estimator is essentially the same as in Q except 
that here it is the full-sky version and the noise Nlm is 
not obtained from the phase-randomized data. We sub- 
tract the Gaussian bias for each realization of the estima- 
tor, and all the estimated power spectra are averaged to 
get the binned power spectra {C£ d ) for the 6-th bin [lfij]. 
The averaged power spectrum in a range of L labeled by 
the index b is 

ndd _ LyL + 1) dd 

C » -Z. b{b+1) °l ■ ( 29 ) 
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FIG. 3: WMAP noise for each DA and the TT power spec- FIG. 5: The averaged reconstruction including noise (C£ st ) 
trum as a function of L. (blue) of WMAP data and the Gaussian bias Nf ] (red) from 

700 realizations. Since lensing is approximately 100 times 
smaller than C£ 3t , the two curves are almost indistinguish- 
able; however, this confirms the precision of the noise model. 
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FIG. 4: Comparison of Al (Eq. l|27[)) and the expected lens- 
ing signal as function of L. The estimator noise is about two 
orders of magnitude higher that the signal Cjf , indicating the 
difficulty of detecting lensing from WMAP-7 data alone. 



FIG. 6: The averaged reconstruction including noise (C£ st ) 
(blue) of simulated WMAP data and the Gaussian bias 
(red) from 700 realizations. Since lensing is approximately 
100 times smaller than C£ st , the two curves are almost indis- 
tinguishable; however, this confirms the precision of the noise 
model. 
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FIG. 7: The reconstructed power spectra (Cf d ) of the de- 
flection angle field from all correlations of WMAP's W- and 
V-band DAs. The green curve is the simulated lensing signal, 
and the data points are the reconstructed lensing signal from 
simulations (red), and the reconstructed lensing signal from 
data (blue). The red and blue data points show the consis- 
tency between the simulated and real WMAP data for the 
lensing reconstruction. 
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FIG. 8: The normalized likelihood distribution for C for all 
21 correlations of WMAP's W- and V-band DAs. 



FIG. 9: Curl null test for all correlations of WMAP's W- and 
V-band DAs: Cf 5 from the simulated WMAP data (red), 
and C S L S from the real WMAP data (blue), for comparison, 
the simulated lensing signal Cf d (solid green). The red and 
blue data points show the consistency between the simulated 
and the real WMAP data for the curl null test. 



The statistical uncertainty is given by <jb = [((C& — 
Cb) 2 )] 5 • After the subtraction of the Gaussian bias, there 
remains the higher order biases, see 0] (where it was 
called "null bias"), 0, and 0. 
We expand T in harmonic space as 



Tlm = Tlm + STlm + S Tlm + 6 Tlm 



(30) 



see [2l( . Here the power n in S n denotes the order in 
We expand the noise bias as 



N L = < } 
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(31) 



where the index n in A/w denotes the order of its de- 
pendence upon [</> 2 ]™, excluding terms that contribute 
to the lensed power spectrum. The four-point func- 
tion (d^jvf <t1ljv/) contains terms of different order in 
S n T. A term of the type (STSTTT) contributes to 
C[ d and the first order noise while terms of 

the type (8TSTST5T), (6 2 TS 2 Tff), {S 2 TST8TT}, and 

(5 3 T5TTT) generate the second order noise N] 2 ' . Fol- 
lowing the higher order bias term is calculated as 
the difference between the estimated power spectrum and 
the sum of its prediction and the lowest order noise (i.e., 
Gaussian bias): Cf* - (Cf + N { ° } ), using Monte Carlo 
simulations. 

We study the statistical significance of the detection as 
follows. Following the reconstructed power spectra 
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(jiphs) arc com p arec i w ith their theoretical prior C*( th ) by 
minimizing a % 2 defined as 
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and varying C. Here ^4 or _B label the range in L, and 
Ca or Cb is the band-power. The covariance matrix C 
is calculated from the Monte Carlo simulation as Cab = 



((c ( sim) 
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The best fit C is 



obtained by setting the derivative of \ to zero: 
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A non-zero value of C indicates the presence of lensing. 
The signficance of a non-zero value can be judged if its 
variance is known. The variance of C is given by 



(AC) 2 = 
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and the significance of the detection of lensing is C/AC. 

We show the higher order bias in Figure [TJ The higher 
order bias N L ^ ± N L 2 ^ ± ... is seen to be negative for 
L < 20 and positive for L > 170 and consistent with zero 
for 20 < L < 170 where the amplitude is -0.42 ± 0.98 
(0.43er), compared to the simulated lensing signal Cff by 
using 15 bins with AL = 10 starting from L = 20. In 
Figure [21 the likelihood of the amplitude of the higher 
order bias limited to the region 20 < L < 170 confirms 
that the bias is consistent with zero. Thus subtraction of 
the higher order bias is not required as long as we limit 
L to this region. 



TABLE I: Measurements of lensing C and its significance 
C/AC. 



Tta+'s col 


c 


C 1 AC 


WMAP-7 ALL" 


1.27 ±0.98 


1.30cr 


WMAP-7 V+W 


0.97 ±0.47 


2.06(t 


WMAP-1 ALLxLRGs c 


1.0 ± 1.1 


0.91a 


WMAP-3 ALLx(LRGs+QSOs+NVSS) <i 


1.06 ±0.42 


2.52a 


WMAP-3 (Q±V±W)xNVSS E 


1.15 ±0.34 


3.38a 


ACT' 


1.16 ±0.29 


4.00cr 


SPT 9 




~ 4.90a 



"All 21 correlations of WMAP-7's W- and V-band DAs in this 
work. 

6 [1 WMAP-7 V and W bands. 

c I WMAP-1 W- and V-band DAs, LRGs. 

d I WMAP-3 W- and V-band DAs, LRGs, QSOs and NVSS. 

e 5 1 WMAP-3 Q-, V-, W-band DAs, NVSS. 

f £ | ACT temperature maps. 

9 [9( SPT temperature maps. 



TABLE II: Summary of C and its significance C/AC for this 
work. 



Type 


C 


C/AC 


higher order bias 


-0.42 ±0.98 


0.43a 


curl null test 


0.38 ±0.79 


0.47a 


reconstructed lensing 


1.27 ±0.98 


1.30(7 



magnitude higher than the signal Cjf; indicative of the 
difficulty of extracting the lensing from the noisy data. 
We calculate the noise in each band from WMAP's data 
instead of using an analytical form as [3| and @ do. The 
noise is simulated according to the prescription in [l9j |. 
and the beam transfer functions are supplied by WMAP. 



V. WMAP 7- YEAR DATA 

The lensing reconstruction depends most sensitively on 
the high-i modes which are supplied by WMAP's DAs 
in the V (2 DAs) and W (4 DAs) frequency bands. Thus 
we use WMAP's coadded temperature maps with r9 res- 
olution (Hcalpix's n s ide = 512) using all possible distinct 
pairings: three auto-correlations for the two V-band DAs, 
ten auto-correlations for the four W-band DAs, and eight 
cross-correlations between the W- and V-band DAs for 
a total of 21 correlations (labeled "ALL"). Smith et al. 
@, used the Q-band DAs in addition to the W- and V- 
band DAs of WMAP 3-year temperature maps. Hirata 
ct al. 0, used 153 one- year DAs from the WMAP 3-year 
data in the W- and V-bands. Recently, Smidt et al. Q 
used the W- and V-frequcncy bands of the WMAP 7- 
year data. This work adopts six DAs of WMAP's 7-year 
temperature map, making the data selection slightly dif- 
ferent from other work, although the same signal-to-noise 
is expected. The WMAP temperature maps contain very 
high levels of noise as shown in Figure [31 The normaliza- 
tion factor Al shown in Figure 2] is about two orders of 



VI. SIMULATION AND ANALYSIS 

We use the CAMB code [22| to obtain the 
power spectra Cf T , and Cf using a six parame- 
ter ACDM model with P = ojbh 2 ,u> c h 2 , h,r, A s ,n s = 
0.0226,0.112,0.70,0.09,2.1 x 10~ 9 ,0.96. These are input 
into a pipeline that has elements as follows. 

Gaussian maps of the deflection angle field d(n) and 
unlensed temperature T(n) are generated using their re- 
spective power spectra. We use Cf T ' , and Cf^ to create 
one realization of the simulated deflection field and lensed 
temperature maps T(n) are generated using Eq. ((TJ) with 
the T(n) and d(n) found above. 

Using the inverse SHT, the temperature maps are con- 
verted into harmonic modes a\ m which are convolved 
with the beam transfer functions bi. Using SHT, these 
are transformed into configuration space and WMAP- 
based noise is added. We mask the galactic plane and 
point sources using WMAP's KQ75 mask. Using inverse 
SHT, the resulting maps are transformed back into har- 
monic space where new a/ m are kept up to l max = 750 
and \m max \ = 750. 
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The noise simulation is crucial to this work because the 
(o) 



' is one hundred times larger than Cff . The six V- 
and W-band DAs, labeled by a = VI, V2, Wl, W2, W3, 
WA, have different noise variances, different beam trans- 
fer functions, and different relative phases. To mimic the 
WMAP DAs, we simulate the Gaussian bias as follows. 
Using 

f (8)Q (n) = M(n)[y dn'fW(n')B Q (n,n') 



(35) 



we set the index i (an arbitrary running index) for Eq. 
(|35p and generate an unlensed temperature map (n) , 
and six noise maps A« Q (n), a = VI, V2, Wl, W2, W3, 
WA. Then we make an observed map using Eq. 

(|35|) , and repeat this procedure to make another observed 
map T^(n). Subsequently, we calculate the Gaussian 
bias using Eq. (p?5"|) for the pair (a»/3). In the same 
way, we generate 21 realizations for all the correlations. 
Finally we increase the index i, and repeat the whole 
procedure until the ensemble {N^} has 700 elements. 

We proceed in a similar manner simulating the recon- 
struction including noise, except setting T^(n) = T(n) 
and 7V« Q (n) = N a (n). Using 

T« Q (n) = M(n) J dri T {t \n')B a {n, n') 



+ A^ Q (n) 



(36) 



we set the index i for Eq. Q36p. and generate a lensed 
temperature map T^(n), and six noise maps N^ a (n), 
a = VI, V2, Wl, W2, W3, W4. Then we make an 
observed map using Eq. (j3"6")l and repeat this 

procedure to make another observed map T^'^(n). Sub- 
sequently we calculate the reconstruction including noise 
C£ s * using Eq. (|2~i]) for the pair (a • /?). In the same 
way, we generate 21 realizations for all the correlations. 
Finally we increase the index i, and repeat the whole pro- 
cedure until the ensemble {C£ s *} has 700 elements. Eq. 
is then used to obtain the deflection power spectrum 



We show the reconstruction including noise C£ s * and 

the Gaussian bias in Figures[5j and|6]for the real and 
the simulated WMAP data, respectively. The simulation 
is consistent with the data, and we confirm that the two 
terms in Eq. (|28|) nearly have the same magnitude, and 
the lensing induced difference is not visible because the 
lensing signal Cff is one hundred times smaller than the 
Gaussian bias . We use Eq. (|28| to calculate the 
reconstructed lensing power spectra in Figure [7l 

The likelihood distribution of C is shown in Figure 
[51 where it is seen lensing is detected at only 1.30(7 
confidence level. 



VII. CURL NULL TEST 

To check for systematic effects, we employ the "curl 
null test". The deflection angle field can be written as 
the sum of a gradient and a curl term [23] : 



D 1 (n) = d i (n)+e <3 -V J '<y(n). (37) 
The first term leads to the Hu estimator 0, Q 

A TT 



a LM 
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(38) 

whose efficient form is given in Eq. (|26|) . here o^4 T (n) is 
given by Eq. dHJ and 



C TT 



(39) 



Cf. 



The estimator for the curl part in Eq. (|37|) is 

^ = £e^^L= / dnYl M V t [ A T (n)V ]0 B T (n)] 

(40) 

and the corresponding efficient form is 



?TT{ a .0) 

'LM 



TT(ax/3) 



■\y/L(L+l) 



(3 La / dn +1 Y£ M0 A T ^X^ _ aw / dn^Y^ A T ^Y^ 



A 



TT(/3xa) 



f3 LQ J dn +1 Yl M0 A T ^X^-a LQ J da- X Y£ M A T ^Y^ 



(41) 



which can be compared with Eq. (|26p . We show the re- from the real and the simulated WMAP data separately 
suiting power spectra Cff , averaged from 700 realizations in Figure [H The averaged curl component amplitude is 
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FIG. 10: The convergence behavior. The values of mean am- 
plitude C (red), the error AC (blue), and the detection signif- 
icance C/AC (green) of the reconstructed lensing signal C L 
are plotted for every 10 realizations. It is seen that conver- 
gence is reached after 700 realizations. 



0.38 ± 0.79 consistent with zero as expected, compared 
to the simulated Cf d . 



VIII. RESULTS AND DISCUSSION 

In this work, we have applied the optimal quadratic 
estimator to WMAP-7 temperature maps alone for the 
first time. 

We have monitored the convergence behavior for the 
mean value C, the error AC, and the detection signifi- 
cance C/AC of the reconstructed lensing signal Cjf. We 
find that all these quantities converge after producing 
700 realizations of the reconstructed lensing signal, see 
Figure [TO] We determine the significance of the lens- 
ing detection and find C = 1.27 ± 0.98 (1.30er), while 
Smidt et al. found C = 0.97 ±0.47 (2.06ct). The result is 
shown in Table Q] as well as a comparison with 
our results have been corrected by the sky fraction. We 
find evidence for lensing only at 1.30er, using all correla- 
tions of WMAP-7's W- and V-band DAs. The resulting 
constraint on the lensing amplitude differs from Q and 
this can be explicated from several aspects. In terms of 
the estimator, we use the optimal estimator derived from 
minimum variance principle , rather than the kurtosis 
estimator in Q. We adopt the individual beam transfer 
function associated with each DA, not the averaged one 
for each frequency. We have taken into account the im- 



pact of the higher order bias, afterwards restricting the 
reconstruction in a proper multiple range that marginally 
overlaps with . In terms of the noise model, we esti- 
mate the noise in a way which mimics WMAP, not simply 
generating random underlying skies and associated noises 
with independent phases. All these factors may jointly 
contribute to the difference between us and Smidt et al. 
A summary of various tests in this work is shown in Table 
UT1 We do not observe a significant lensing signal from 
the WMAP 7-year temperature data. 

We did not apply a correction for higher order bias 
terms Nj/, Nj? , because they are expected to be 
small owing to the fact that we limited the region of L to 
20 < L < 170, where the higher order bias is consistent 
with zero. The higher order bias can be obtained via an 
iterative solution [18| but it is computationally demand- 
ing and not warranted in the present case because we do 
not obtain a significant signal. 

We applied the curl null test to all the correlations of 
W- and V-band DAs as a systematic check, since we ob- 
serve a small amount of power from the reconstructed 
gravitational lensing signal (Figure [7]) . The reconstruc- 
tion procedure passes the curl null test. 

The effects of beam systematics and the galactic and 
foreground contaminations are quite small compared to 
the statistical error. We do not correct the statistical 
result for the presence of point sources because they in- 
troduce negligible systematics fl9| . 

We have demonstrated, using a nearly optimal esti- 
mator, that WMAP-7 data does not have the power to 
detect gravitational lensing, which is unfortunate since 
WMAP data is the only publicly available data set with 
sufficient angular resolution to detect lensing. However, 
WMAP-7 does have value as a publicly available tool to 
assess the efficacy of lensing algorithms and to test for 
systematic biases. 
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Gravitational lensing by large scale structure introduces non-Gaussianity into the Cosmic Mi- 
crowave Background and imprints a new observable, which can be used as a cosmological probe. 
We apply a four-point estimator to the Wilkinson Microwave Anisotropy Probe (WMAP) 7-year 
coadded temperature maps alone to reconstruct the gravitational lensing signal. The Gaussian bias 
is simulated and subtracted, and the higher order bias is investigated. We measure a gravitational 
lensing signal with a statistical amplitude of C = 1.27 ± 0.98 using all the correlations of the W- 
and V-band Differencing Assemblies (DAs). We therefore conclude that WMAP 7-year data alone, 
can not detect lensing. 



I. INTRODUCTION 

Gravitational lensing of the Cosmic Microwave Back- 
ground (CMB) provides information on the mass distri- 
bution between the surface of last scattering and the ob- 
server, thus potentially providing information, for exam- 
ple, on dark energy and neutrino masses. In addition, 
gravitational lensing causes E'-modes to be converted into 
large angular scale i?-modes, thereby potentially contam- 
inating B-mode signature of inflationary gravitational 
waves Because lensing deflects CMB photons by ap- 
proximately 3', a pcrturbative treatment to first order is 
generally valid. An estimator for the deflection angle has 
been devised by Hu @, H| . 

The first attempt to detect lensing by Hirata et al. 
[H used the cross-correlation between the WMAP 1-year 
data and selected luminous red galaxies (LRGs) from the 
Sloan Digital Sky Survey (SDSS). No statistically signif- 
icant signal was found. The first detection of lensing 
was performed by Smith et al. Q who used the cross- 
correlation between the NRAO VLA Sky Survey (NVSS) 
of radio galaxies with a higher mean redshift than the 
Sloan LRGs and a fully-optimal lensing estimator on the 
statistically more powerful WMAP 3-year data. Evi- 
dence for lensing was found at the 3Aa level. Using a 
similar estimator as in Hirata et al. @ obtained re- 
sults consistent with, though at slightly lower significance 
than Q, using WMAP 3-year data, LRGs and quasars 
from the SDSS data, as well as data from the NVSS. Rc- 
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cently, Smidt et al. Q used an estimator based upon 
the kurtosis of the CMB temperature four-point corre- 
lation function to estimate lensing from WMAP 7-year 
data only and claimed evidence for lensing at the 2a level. 
Recently, the Atacama Cosmology Telescope (ACT) col- 
laboration successfully detected gravitational lensing Q 
at the 4cr level. The South Pole Telescope (SPT) de- 
tected the effects of gravitational lensing on the angular 
power spectrum 0. 

In this paper we present a search for gravitational lens- 
ing using the WMAP 7-year data alone and the standard 
optimal quadratic estimator 0, Q which differs from the 
kurtosis estimator of Q ■ We apply the quadratic estima- 
tor to WMAP-7 temperature maps alone for the first time 
in the hopes that our analysis might serve as a touchstone 
allowing for consistent comparison between different lens- 
ing extraction techniques. We review the notation for 
full-sky reconstruction of gravitational lensing in Section 

II. We discuss the sky-cut used in our analysis in Section 

III. Then we introduce our modified estimator in Sec- 
tion IV making use of the optimal quadratic estimator 
of 0. We introduce the WMAP 7 -year data in Section 
V, and describe the details of the calculations, including 
the noise model, and analysis in Section VI. Results of 
a null test are shown in Section VII, and we discuss the 
conclusions of our work in Section VIII. 



II. GRAVITATIONAL LENSING 

The effect of lensing on the CMB's primordial temper- 
ature T in direction n can be represented by 



T(n) = f(n + d(n)), 



(1) 



2 



where T is the lensed temperature and d(n) = V</>, with 
4> being the lensing potential. The two-point correlation 
function of the temperature field following is: 
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where the second term encodes the effects of lensing with 
the weighting factor f^J, given by 

fiui' = Cf T oFl'U + Cf T a Fi L i>- (3) 

Here Cj T are the unlensed temperature power spectra, 
and 



qFiw = 



(2/ + l)(2/ / + l)(2£+l) 

47T 



l[L(L+l) + l'(l' + l)-l(l + l)}( l Q J [J) 



(4) 



The lensing estimator is constructed from an average 
over a pair of two-point correlations 0, Q and has the 
form 
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The requirement that the estimator in Eq. 
biased and has minimal variance results in 
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L(L + 1)(2L + 1) 



and 
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(6) 



(7) 



with C\ ot = Cf T +Nf T , where Cf T are the lensed power 
spectra and Nf T is the instrumental noise. In the fol- 
lowing, the summations are from / and V = to 750 
and |m| < i, \m'\ < I'. The WMAP 7-year data do not 
contain additional information at higher multipoles. 

To reduce computation time we follow (Tol | and define 
three maps for the TT estimator: 



A T (n)=Y,-^tT lrn0 Y lm (n), 



I in 



(8) 



and take the inverse Spherical Harmonic Transform 
(SHT) of A T X and A T Y to get 



Tlm = Plo I dn +1 Y£ M A T X 



with 



-(2) 
- LM 



ft, 
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(11) 



(12) 
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(14) 



Using Eqs. (jHJl, © and (fl0|) the expression for d^Jt in 
Eq. ([5]) becomes 



JTT _ 
a LM — 



r T (l) T (2) 1 



(15) 



A similar procedure is followed for the efficient calcula- 
tion of Ajr in Eq. (j6|) . The resulting expression is given 
in [ll| (originally proposed in fl2l|): 
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with the £ T given by 
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ff ±1 (g) = g + 1 ) C i?+ ^TT ^lW. (19) 



here d' ss / (0) are Wigner d-functions. 



III. SKY CUT 



C TT 

X ( n ) = E -^t T ^ a io +iY lm (n), 

Irn ' 



F ( n )=E ^SF^mAo-l^m(n), 



In order to eliminate contaminated data, regions such 
(9) as the galactic plane and bright point sources in the full- 
sky map must be removed using a mask, thereby intro- 
ducing a sky-cut. For example, in 0, the Kp2 mask 
was used to make 84.7% of the sky uncontaminated. In 
(10) 0; the more conservative KQ75 mask was used to clean 
artifacts around the galactic plane and point sources. 
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The sky-cut can be removed as a separate component 
to get a full-sky map before we process the data. One 
such technique is the "inpainting" method in which the 
estimated values of pixels in the map are substituted for 
those removed by the mask. Perotto et al. have simu- 
lated the full sky reconstruction for PLANCK Pj. The 
full-sky map recovered in this way will bias the lensing 
reconstruction slightly. 

Another method proposed by A. Benoit-Levy [l4[ 
apodizes the masked regions of the map and inpaints the 
masked regions of the map by constrained Gaussian ran- 
dom values of the unlensed temperature. In this way, 
the sky-cut-induced coupling approximately reduces to a 
unit matrix. However, for WMAP, we have to remove a 
big portion of the sky, reducing / s k y dramatically to 0.3. 
The unbiased estimator could be scaled up by a factor 
of l// s ky, but the signal-to-noise ratio would be reduced 
significantly. This means the uncertainty of the recon- 
structed signal would be larger. 

As opposed to a separate-component solution, we ob- 
tain an all-inclusive lensing reconstruction pipeline, using 
the built-in filter of the estimator to treat the data with- 
out pre-conditioning it. The optimal estimator for the 
potential based on the maximum likelihood is derived by 
Hirata fl5j . The full inverse variance (C + N)- 1 , instead 
of (Cj T + Nj / T )~ 1 , was used by @ because it is an op- 
timal filter when there are sky-cuts and inhomogeneous 
noise. The sky-cut generates artifacts in harmonic space, 
as does lensing. (C + N) _1 can be used to filter those 
modes affected by the sky-cut. However, we do not use 
this filter because the inversion of (C + N) is computa- 
tionally challenging [f| ; instead we use the estimator Eq. 
([5]) which is identical to the one of [f| , and it is an excel- 
lent approximation to the maximum likelihood estimator. 
We note that, while (Cj T + A r J T ) _1 will be suboptimal 
to a full (C + N) -1 filter, it preserves the simplicity and 
efficiency of the lensing reconstruction procedure. 
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FIG. 1: The higher order bias calculated from {CT - iv£ 0) ) - 
Ct d for all correlations of the WMAP's W- and V-band DAs. 
The simulated higher order bias from averaging 700 (to be 
discussed in Figure I10p realizations is shown in orange. For 
comparison, the simulated lensing signal is shown in green. 
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IV. THE LENSING ESTIMATOR 



following Eq. , and a harmonic mode of the Gaussian 
bias is estimated as 



For WMAP, we modify the estimator slightly to deal 
with the instruments' anisotropic temperature noise. 
The observed lensed temperature map T is given by 

T(n) = M(n) [ J dn'T(n')B(n, n') + 7V(n)] (20) 

and likewise the "observed" unlensed temperature map 
f is 

T(n) = M(n) [ J drif(n')B(n, n') + N(n)] (21) 

Here M(n) represents the mask, i?(n,n') the beam, and 
N(n) the noise. 

For a pair of maps a and f3, "TT(a x /?)" denotes the 
cross-correlation between these two temperature maps. 
A harmonic mode of the reconstruction including noise 
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Here C are the power spectra of the observed lensed 
temperature, determined from (T; TO T;' m '). As was done 
in [8( and [16j we use the same power spectra in Eq. 
([22]) and Eq. ([23"]). In order to deal with the non- 
uniform noise distribution in the WMAP data, we sym- 
metrize d^^} QX ' 9 '' as in 0, denoting the symmetrized 
cross-correlation u TT(a» (3)" between these two temper- 
ature maps, 
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and including noise, and = {N^ m Nlm) as the Gaussian 

, „. rprp.H n bias, with the superscript "TT(a»/3)" omitted. Thus we 

lm +jn lm ^5) obtain 

We refer to C£ s * = (d* LM d,LM) a s the reconstruction 
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FIG. 2: The normalized likehood of the amplitude of the 
higher order bias limited to the region 20 < L < 170, to 
the simulated lensing signal. This confirms that the higher 
order bias is consistent with zero and negligible. 



following a reasoning similar to the one near the end of 
Section II. 

The two-point correlation of the Gaussian bias esti- 
mator is essentially a four-point correlation function of 
the primordial temperature modes. It should be care- 
fully subtracted since, for a noise-dominated experiment 
such as WMAP, the Gaussian four-point bias is sev- 



eral orders of magnitude larger than the lensing power 
spectra. In [§[ phase-randomized data maps are used 
to simulate this Gaussian bias. However, this approach 
does not work for the present lensing reconstruction since 
WMAP's noise is not isotropic. Evidence for this can be 
seen from the normalization factor A^ T< ' ax ^ which is 
not equal to j\f^ TT ^ a '^ whereas they should be equal 
for isotropic noise (Io| . The normalization factor Eq. 
(|27p only contains the partial contribution coming from 
the non-isotropic noise while the Gaussian bias squared 
from Eq . (|2~5j) consists of all the correlations generated by 
the non-isotropic noise, see [T3| and [l8[ . If the phases of 
the WMAP temperature maps are randomized in order 
to remove the lensing-induccd coupling between modes, 
it will also remove the strong correlation of the noise. 
The Gaussian bias calculated in this way will be signifi- 
cantly lower than that from the standard approach [l9j |. 
So we have to perform simulations which use the simu- 
lated WMAP noise and temperature maps, rather than 
the randomized WMAP data to get the Gaussian bias 
term. 

The deflection power spectrum is 

ct = < [ d £<~»] * - [<T" ] ] * Kit^y 

n (28) 

This estimator is essentially the same as in Q except 
that here it is the full-sky version and the noise Nlm is 
not obtained from the phase-randomized data. We sub- 
tract the Gaussian bias for each realization of the estima- 
tor, and all the estimated power spectra are averaged to 
get the binned power spectra {C£ d ) for the 6-th bin [lfij]. 
The averaged power spectrum in a range of L labeled by 
the index b is 

ndd _ LyL + 1) dd 

C » -Z. b{b+1) °l ■ ( 29 ) 
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FIG. 3: WMAP noise for each DA and the TT power spec- FIG. 5: The averaged reconstruction including noise (C£ st ) 
trum as a function of L. (blue) of WMAP data and the Gaussian bias Nf ] (red) from 

700 realizations. Since lensing is approximately 100 times 
smaller than C£ 3t , the two curves are almost indistinguish- 
able; however, this confirms the precision of the noise model. 
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FIG. 4: Comparison of Al (Eq. l|27[)) and the expected lens- 
ing signal as function of L. The estimator noise is about two 
orders of magnitude higher that the signal Cjf , indicating the 
difficulty of detecting lensing from WMAP-7 data alone. 



FIG. 6: The averaged reconstruction including noise (C£ st ) 
(blue) of simulated WMAP data and the Gaussian bias 
(red) from 700 realizations. Since lensing is approximately 
100 times smaller than C£ st , the two curves are almost indis- 
tinguishable; however, this confirms the precision of the noise 
model. 
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FIG. 7: The reconstructed power spectra (Cf d ) of the de- 
flection angle field from all correlations of WMAP's W- and 
V-band DAs. The green curve is the simulated lensing signal, 
and the data points are the reconstructed lensing signal from 
simulations (red), and the reconstructed lensing signal from 
data (blue). The red and blue data points show the consis- 
tency between the simulated and real WMAP data for the 
lensing reconstruction. 



0.5 




FIG. 8: The normalized likelihood distribution for C for all 
21 correlations of WMAP's W- and V-band DAs. 



FIG. 9: Curl null test for all correlations of WMAP's W- and 
V-band DAs: Cf 5 from the simulated WMAP data (red), 
and C S L S from the real WMAP data (blue), for comparison, 
the simulated lensing signal Cf d (solid green). The red and 
blue data points show the consistency between the simulated 
and the real WMAP data for the curl null test. 



The statistical uncertainty is given by <jb = [((C& — 
Cb) 2 )] 5 • After the subtraction of the Gaussian bias, there 
remains the higher order biases, see 0] (where it was 
called "null bias"), 0, and 0. 
We expand T in harmonic space as 



Tlm = Tlm + STlm + S Tlm + 6 Tlm 



(30) 



see [2l( . Here the power n in S n denotes the order in 
We expand the noise bias as 



N L = < } 



N 



(i) 



N 



(2) 



(31) 



where the index n in A/w denotes the order of its de- 
pendence upon [</> 2 ]™, excluding terms that contribute 
to the lensed power spectrum. The four-point func- 
tion (d^jvf <t1ljv/) contains terms of different order in 
S n T. A term of the type (STSTTT) contributes to 
C[ d and the first order noise while terms of 

the type (8TSTST5T), (6 2 TS 2 Tff), {S 2 TST8TT}, and 

(5 3 T5TTT) generate the second order noise N] 2 ' . Fol- 
lowing the higher order bias term is calculated as 
the difference between the estimated power spectrum and 
the sum of its prediction and the lowest order noise (i.e., 
Gaussian bias): Cf* - (Cf + N { ° } ), using Monte Carlo 
simulations. 

We study the statistical significance of the detection as 
follows. Following the reconstructed power spectra 
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(jiphs) arc com p arec i w ith their theoretical prior C*( th ) by 
minimizing a % 2 defined as 
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and varying C. Here ^4 or _B label the range in L, and 
Ca or Cb is the band-power. The covariance matrix C 
is calculated from the Monte Carlo simulation as Cab = 



((c ( sim) 
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The best fit C is 



obtained by setting the derivative of \ to zero: 
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A non-zero value of C indicates the presence of lensing. 
The signficance of a non-zero value can be judged if its 
variance is known. The variance of C is given by 



(AC) 2 = 
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and the significance of the detection of lensing is C/AC. 

We show the higher order bias in Figure [TJ The higher 
order bias N L ^ ± N L 2 ^ ± ... is seen to be negative for 
L < 20 and positive for L > 170 and consistent with zero 
for 20 < L < 170 where the amplitude is -0.42 ± 0.98 
(0.43er), compared to the simulated lensing signal Cff by 
using 15 bins with AL = 10 starting from L = 20. In 
Figure [21 the likelihood of the amplitude of the higher 
order bias limited to the region 20 < L < 170 confirms 
that the bias is consistent with zero. Thus subtraction of 
the higher order bias is not required as long as we limit 
L to this region. 



TABLE I: Measurements of lensing C and its significance 
C/AC. 



Tta+'s col 


c 


C 1 AC 


WMAP-7 ALL" 


1.27 ±0.98 


1.30cr 


WMAP-7 V+W 


0.97 ±0.47 


2.06(t 


WMAP-1 ALLxLRGs c 


1.0 ± 1.1 


0.91a 


WMAP-3 ALLx(LRGs+QSOs+NVSS) <i 


1.06 ±0.42 


2.52a 


WMAP-3 (Q±V±W)xNVSS E 


1.15 ±0.34 


3.38a 


ACT' 


1.16 ±0.29 


4.00cr 


SPT 9 




~ 4.90a 



"All 21 correlations of WMAP-7's W- and V-band DAs in this 
work. 

6 [1 WMAP-7 V and W bands. 

c I WMAP-1 W- and V-band DAs, LRGs. 

d I WMAP-3 W- and V-band DAs, LRGs, QSOs and NVSS. 

e 5 1 WMAP-3 Q-, V-, W-band DAs, NVSS. 

f £ | ACT temperature maps. 

9 [9( SPT temperature maps. 



TABLE II: Summary of C and its significance C/AC for this 
work. 



Type 


C 


C/AC 


higher order bias 


-0.42 ±0.98 


0.43a 


curl null test 


0.38 ±0.79 


0.47a 


reconstructed lensing 


1.27 ±0.98 


1.30(7 



magnitude higher than the signal Cjf; indicative of the 
difficulty of extracting the lensing from the noisy data. 
We calculate the noise in each band from WMAP's data 
instead of using an analytical form as [3| and @ do. The 
noise is simulated according to the prescription in [l9j |. 
and the beam transfer functions are supplied by WMAP. 



V. WMAP 7- YEAR DATA 

The lensing reconstruction depends most sensitively on 
the high-i modes which are supplied by WMAP's DAs 
in the V (2 DAs) and W (4 DAs) frequency bands. Thus 
we use WMAP's coadded temperature maps with r9 res- 
olution (Hcalpix's n s ide = 512) using all possible distinct 
pairings: three auto-correlations for the two V-band DAs, 
ten auto-correlations for the four W-band DAs, and eight 
cross-correlations between the W- and V-band DAs for 
a total of 21 correlations (labeled "ALL"). Smith et al. 
@, used the Q-band DAs in addition to the W- and V- 
band DAs of WMAP 3-year temperature maps. Hirata 
ct al. 0, used 153 one- year DAs from the WMAP 3-year 
data in the W- and V-bands. Recently, Smidt et al. Q 
used the W- and V-frequcncy bands of the WMAP 7- 
year data. This work adopts six DAs of WMAP's 7-year 
temperature map, making the data selection slightly dif- 
ferent from other work, although the same signal-to-noise 
is expected. The WMAP temperature maps contain very 
high levels of noise as shown in Figure [31 The normaliza- 
tion factor Al shown in Figure 2] is about two orders of 



VI. SIMULATION AND ANALYSIS 

We use the CAMB code [22| to obtain the 
power spectra Cf T , and Cf using a six parame- 
ter ACDM model with P = ojbh 2 ,u> c h 2 , h,r, A s ,n s = 
0.0226,0.112,0.70,0.09,2.1 x 10~ 9 ,0.96. These are input 
into a pipeline that has elements as follows. 

Gaussian maps of the deflection angle field d(n) and 
unlensed temperature T(n) are generated using their re- 
spective power spectra. We use Cf T ' , and Cf^ to create 
one realization of the simulated deflection field and lensed 
temperature maps T(n) are generated using Eq. ((TJ) with 
the T(n) and d(n) found above. 

Using the inverse SHT, the temperature maps are con- 
verted into harmonic modes a\ m which are convolved 
with the beam transfer functions bi. Using SHT, these 
are transformed into configuration space and WMAP- 
based noise is added. We mask the galactic plane and 
point sources using WMAP's KQ75 mask. Using inverse 
SHT, the resulting maps are transformed back into har- 
monic space where new a/ m are kept up to l max = 750 
and \m max \ = 750. 



8 



The noise simulation is crucial to this work because the 
(o) 



' is one hundred times larger than Cff . The six V- 
and W-band DAs, labeled by a = VI, V2, Wl, W2, W3, 
WA, have different noise variances, different beam trans- 
fer functions, and different relative phases. To mimic the 
WMAP DAs, we simulate the Gaussian bias as follows. 
Using 

f (8)Q (n) = M(n)[y dn'fW(n')B Q (n,n') 



(35) 



we set the index i (an arbitrary running index) for Eq. 
(|35p and generate an unlensed temperature map (n) , 
and six noise maps A« Q (n), a = VI, V2, Wl, W2, W3, 
WA. Then we make an observed map using Eq. 

(|35|) , and repeat this procedure to make another observed 
map T^(n). Subsequently, we calculate the Gaussian 
bias using Eq. (p?5"|) for the pair (a»/3). In the same 
way, we generate 21 realizations for all the correlations. 
Finally we increase the index i, and repeat the whole 
procedure until the ensemble {N^} has 700 elements. 

We proceed in a similar manner simulating the recon- 
struction including noise, except setting T^(n) = T(n) 
and 7V« Q (n) = N a (n). Using 

T« Q (n) = M(n) J dri T {t \n')B a {n, n') 



+ A^ Q (n) 



(36) 



we set the index i for Eq. Q36p. and generate a lensed 
temperature map T^(n), and six noise maps N^ a (n), 
a = VI, V2, Wl, W2, W3, W4. Then we make an 
observed map using Eq. (j3"6")l and repeat this 

procedure to make another observed map T^'^(n). Sub- 
sequently we calculate the reconstruction including noise 
C£ s * using Eq. (|2~i]) for the pair (a • /?). In the same 
way, we generate 21 realizations for all the correlations. 
Finally we increase the index i, and repeat the whole pro- 
cedure until the ensemble {C£ s *} has 700 elements. Eq. 
is then used to obtain the deflection power spectrum 



We show the reconstruction including noise C£ s * and 

the Gaussian bias in Figures[5j and|6]for the real and 
the simulated WMAP data, respectively. The simulation 
is consistent with the data, and we confirm that the two 
terms in Eq. (|28|) nearly have the same magnitude, and 
the lensing induced difference is not visible because the 
lensing signal Cff is one hundred times smaller than the 
Gaussian bias . We use Eq. (|28| to calculate the 
reconstructed lensing power spectra in Figure [7l 

The likelihood distribution of C is shown in Figure 
[51 where it is seen lensing is detected at only 1.30(7 
confidence level. 



VII. CURL NULL TEST 

To check for systematic effects, we employ the "curl 
null test". The deflection angle field can be written as 
the sum of a gradient and a curl term [23] : 



D 1 (n) = d i (n)+e <3 -V J '<y(n). (37) 
The first term leads to the Hu estimator 0, Q 

A TT 



a LM 



y/L(L + l) 



dnYl M ^[ A T (n)V lQ B T (n)] 



(38) 

whose efficient form is given in Eq. (|26|) . here o^4 T (n) is 
given by Eq. dHJ and 



C TT 



(39) 



Cf. 



The estimator for the curl part in Eq. (|37|) is 

^ = £e^^L= / dnYl M V t [ A T (n)V ]0 B T (n)] 

(40) 

and the corresponding efficient form is 



?TT{ a .0) 

'LM 



TT(ax/3) 



■\y/L(L+l) 



(3 La / dn +1 Y£ M0 A T ^X^ _ aw / dn^Y^ A T ^Y^ 



A 



TT(/3xa) 



f3 LQ J dn +1 Yl M0 A T ^X^-a LQ J da- X Y£ M A T ^Y^ 



(41) 



which can be compared with Eq. (|26p . We show the re- from the real and the simulated WMAP data separately 
suiting power spectra Cff , averaged from 700 realizations in Figure [H The averaged curl component amplitude is 
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T 




O.o I — ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 

100 200 300 400 500 600 700 

number of iteration 



FIG. 10: The convergence behavior. The values of mean am- 
plitude C (red), the error AC (blue), and the detection signif- 
icance C/AC (green) of the reconstructed lensing signal C L 
are plotted for every 10 realizations. It is seen that conver- 
gence is reached after 700 realizations. 



0.38 ± 0.79 consistent with zero as expected, compared 
to the simulated Cf d . 



VIII. RESULTS AND DISCUSSION 

In this work, we have applied the optimal quadratic 
estimator to WMAP-7 temperature maps alone for the 
first time. 

We have monitored the convergence behavior for the 
mean value C, the error AC, and the detection signifi- 
cance C/AC of the reconstructed lensing signal Cjf. We 
find that all these quantities converge after producing 
700 realizations of the reconstructed lensing signal, see 
Figure [TO] We determine the significance of the lens- 
ing detection and find C = 1.27 ± 0.98 (1.30er), while 
Smidt et al. found C = 0.97 ±0.47 (2.06ct). The result is 
shown in Table Q] as well as a comparison with 
our results have been corrected by the sky fraction. We 
find evidence for lensing only at 1.30er, using all correla- 
tions of WMAP-7's W- and V-band DAs. The resulting 
constraint on the lensing amplitude differs from Q and 
this can be explicated from several aspects. In terms of 
the estimator, we use the optimal estimator derived from 
minimum variance principle , rather than the kurtosis 
estimator in Q. We adopt the individual beam transfer 
function associated with each DA, not the averaged one 
for each frequency. We have taken into account the im- 



pact of the higher order bias, afterwards restricting the 
reconstruction in a proper multiple range that marginally 
overlaps with . In terms of the noise model, we esti- 
mate the noise in a way which mimics WMAP, not simply 
generating random underlying skies and associated noises 
with independent phases. All these factors may jointly 
contribute to the difference between us and Smidt et al. 
A summary of various tests in this work is shown in Table 
UT1 We do not observe a significant lensing signal from 
the WMAP 7-year temperature data. 

We did not apply a correction for higher order bias 
terms Nj/, Nj? , because they are expected to be 
small owing to the fact that we limited the region of L to 
20 < L < 170, where the higher order bias is consistent 
with zero. The higher order bias can be obtained via an 
iterative solution [18| but it is computationally demand- 
ing and not warranted in the present case because we do 
not obtain a significant signal. 

We applied the curl null test to all the correlations of 
W- and V-band DAs as a systematic check, since we ob- 
serve a small amount of power from the reconstructed 
gravitational lensing signal (Figure [7]) . The reconstruc- 
tion procedure passes the curl null test. 

The effects of beam systematics and the galactic and 
foreground contaminations are quite small compared to 
the statistical error. We do not correct the statistical 
result for the presence of point sources because they in- 
troduce negligible systematics fl9| . 

We have demonstrated, using a nearly optimal esti- 
mator, that WMAP-7 data does not have the power to 
detect gravitational lensing, which is unfortunate since 
WMAP data is the only publicly available data set with 
sufficient angular resolution to detect lensing. However, 
WMAP-7 does have value as a publicly available tool to 
assess the efficacy of lensing algorithms and to test for 
systematic biases. 
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